An exactly solvable model of a superconducting 
to rotational phase transition 

D.J. Rowe, C. Bahri, and W. Wijesundera 

Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada 

(February 9, 2008) 



Abstract 

We consider a many-fermion model which exhibits a transition from a super- 
conducting to a rotational phase with variation of a parameter in its Hamilto- 
nian. The model has analytical solutions in its two limits due to the presence 
of dynamical symmetries. However, the symmetries are basically incompatible 
with one another which means that no simple solution exists in intermediate 
situations. Exact (numerical) solutions are nevertheless possible on a com- 
puter and enable one to study the behaviour of competing but incompatible 
symmetries and the phase transitions that result in a semi-realistic situation. 
In spite of the complexity of the mixed-symmetry dynamics, the results are 

remarkably simple and shed light on the nature of phase transitions. 
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Much can be learned about the vahdity of approximate many-body theories by apply- 
ing them to exactly solvable models [jl],^. We consider a model with two simple limits 
(phases): one superconducting; the other rotational. Systems of this kind are prevalent 
in nuclear physics. Indeed, most open-shell nuclei are regarded as being vibrational and 
superconducting or deformed and rotational. Some nuclei appear to be both rotational and 
superconducting. Thus, we consider a model with Hamiltonian 

H = aHi + {l- a)H2 , (1) 

where Hi and H2 are, respectively, Hamiltonians for rotational motion and superconductiv- 
ity; a is a parameter in the range < a < 1. With variation of a, the system undergoes a 
transition from one phase to the other. 

A common approach to understanding such systems is with the Hartree-Fock-Bogolyubov 
approximation or one of its phenomenological variations. Such mean-field theories yield 
many solutions in general: some spherical and superconducting; others deformed and rota- 
tional. Moreover, with variation of the particle number or parameters in the Hamiltonian, 
phase transitions occur in the solution which lies lowest in energy. This motivates one to in- 
quire if the phase transitions are real or an artifact of the HFB constraints. When improved 
approximations are made, e.g., by restoring broken rotational and particle number symme- 
tries by angular momentum and number projection, it appears that the phase transitions 
are smoothed out [^. Remnants may survive and even become sharp in large systems. But, 
one needs to compare the results with exact calculations, to know what they have to say 
about physical phase transitions. 

Exact results for truncated spaces can be obtained with a shell model code. For example, 
Bahri et al. have investigated the competition between superconductivity and rotational 
dynamics for a system of identical nucleons in the nuclear fp shell. The remarkable fact is 
that superconductivity and rotations, separately, can be treated rather well as a result of 
their dynamical symmetries [Q. However, for a combination of the two, one needs a huge 
space to include all the relevant correlations. This is because the dynamical symmetries for 
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rotations and superconductivity are incompatible according to the following definition. 

If a Hilbert space H carries reducible representations of two groups, the actions of the 
two groups are said to be incompatible if they do not commute with one another and if H has 
no proper subspace that is invariant under both groups. For example two Lie group actions 
might be incompatible because there exists no finite-dimensional Lie group that contains 
both groups as subgroups. 

Microscopic models for the rotational states of an A-particle nucleus are given by Hamil- 
tonians of the form 



A A 

H^^ = ^Y.Pn + \rnu' E + ^(^) ' (2) 



n n 



where V{Q) is a low order polynomial in the quadrupole moments of the nucleus The 
building blocks of such a Hamiltonian (^p^, ^ r^, and the components of Q) are elements of 
the (noncompact) symplectic algebra Sp(3,]R). Thus, i^RM is diagonalizable on a computer 
(to within round-off errors). 

Models of superconductivity [§] are given by Hamiltonians of the form 

-f^sc = E ^j<^]m"'jm - E Gjjia^^a^^ajifh'O'j'm' , (3) 

where CLjm} are creation and annihation operators for fermions of angular momentum 
j and 2;-component m; jrh label the time-reverse of the state labelled by jm. Such a 
Hamiltonian has a spectrum generating algebra (SGA) comprising a direct sum of SU(2) 
(quasi-spin) algebras 0. Consequently, it can be diagonalized numerically (if the number 
of single particle states is not too large) [1T0| . 

The actions of the Sp(3, M) and quasispin Lie algebras are incompatible. As a result, it 
is not possible to diagonalize a linear combination of i^RM and Hsc without an unjustified 
truncation of the space. It is then an interesting challenge to find approximate methods for 
determining the spectrum of such combinations. 

To simplify the problem, we restrict the symplectic model to a single spherical har- 



monic oscillator shell. The model then reduces to Elliott's model []IT| for which a suitable 
Hamiltonian is 



H^ = -xQ-Q. (4) 

where Q is the SU(3) quadrupole tensor. This Hamiltonian has spectrum given by 

^(A/iJ) = -xC2(A/i) + 3xJ( J + 1) , (5) 

where (A/i) is an SU(3) highest weight, J is the angular momentum, and 

C2(A/i) =4(A2 + ;u2 + A/i + 3A + 3/i) (6) 

is the value of the SU(3) Casimir invariant. 

A simplification of ifsc is obtained by putting all single-particle energies equal and setting 
Gjji = G. The Hamiltonian H^q then reduces to 

H2 = -GS+S. , (7) 

where 

5- = E,™>o(-l)^'+™«.,-«.-- (8) 

This Hamiltonian has analytical solutions because, as shown by Kerman et al. S± are the 
raising and lowering operators of a single SU(2) algebra with [5*+, SJ\ = 2So and [Sq, S±] = 
±S±. It has a spectrum given by 

E{N, s) = -Gs{s + 1) + lG{N - 2s){N - 2s - 2) . (9) 

With these simplifications, the dynamical symmetries of Hi and H2 can be combined. 
The smallest Lie algebra that contains both SU(2) and SU(3), acting within the space of a 
many-particle spherical harmonic oscillator shell, is the compact symplectic algebra Sp(n); 
an algebra of rank n equal to the number, n = ^(z/ + l)(z/ + 2), of orbital single-particle 
states for the shell with u harmonic oscillator quanta. It is spanned by operators 
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km 

Bij = '^^{ — l)'^^^akimO'kj-m, (10) 
km 

~ 2 ''^^^i'^\im'^kjm — O^k-jmO^kim) ' 



km 

where m is summed over spin values ±1/2, i and j run over the range 1, . . . , z/, and k is 
summed over k fermion types; e.g., k = 2 for neutrons and protons. 

For n > 3, Sp(n) contains EUiott's SU(3) and Kerman's SU(2) quasispin algebra, with 

S+ = \ ^ ^ An , 5*0 = ^ ^ ] C*M , "S*- = ^ ^ ] -Bji , (11) 

i i i 

as subalgebras. Thus, sp(n) is a SGA for the Hamiltonian of eqn. (|l]). We consider here 
the smallest Sp(n) Lie algebra that contains both SU(3) and SU(2), i.e., Sp(3). One recalls 
that 30(5) symmetry (isomorphic to sp{z)) has been used by Zhang and Demler to 
explore the competition between antiferromagnetism and d-wave superconductivity within 
the framework of the Hubbard model. 

The relevant irrep of Sp(3), for present purposes, is the one whose lowest weight state is 
the fermion vacuum state |0). This state satisfies the equations 

5,,|0) = 0, Ci,\Q) = -5,.jk\Q) , (12) 

which means that it is a lowest weight state for an Sp(3) irrep with weight ( 

Basis states for an irrep with this lowest weight are generated by repeated application 
of the Sp(3) raising operators {Aij} to the vacuum state |0). These operators transform as 
a basis for a u(3) irrep of highest weight {200}. Thus, an orthonormal basis {\nKJM)} for 
the desired Sp(3) irrep can be constructed, which reduces the subalgebra chain 

sp{3) D u(3) D su(3) D so(3) D 5o(2) , (13) 

where n = {ni,n2,n^}, a triple of integers with rii > n2 > n^, is a u(3) highest weight 
and JM are the usual so(3) D so(2) angular momentum quantum numbers; K indexes 
multiplicities of J in the SU(3) irrep with highest weight (A,;u) = {rii — ^2,^2 — n^). 



The matrices of such irreps of the symplectic Lie algebras are known |ll3| , p^ ; exphcit 
expressions will be given in a more complete report to follow. 

For K = 2, one obtains Sp(3) representations appropriate for two kinds of nucleon (neu- 
trons and protons) filling the Ip shell. This space is too small; the most deformed SU(3) 
irrep it contains is of highest weight (A/i) = (4, 0) whereas the irreps relevant for the descrip- 
tion of rotational bands in rare-earth nuclei ||15| have A ~ 70 — 100, /i ~ — 10. Such SU(3) 



irreps are found in the nuclear shell model for two kinds of nucleon (i.e., k = 2) filling the 
z/ = 4 and 5 shells. However, mixed SU(3) calculations in these shells are impossible. We 
therefore generate the desired SU(3) irreps in our model, in which the most deformed SU(3) 
irrep is of highest weight {\^) = (2k, 0), by artificially assigning k, values in the 20 — 40 
range. 

We stress that our objective in this paper is not to fit any particular data. It is rather to 
construct an exactly solvable model that can be used to test ideas and many-body theories. 
It should also be understood that our model can be derived in many ways. For example, one 
can regard its Sp(3) algebra as a subalgebra of a larger Sp(n) Lie algebra using pseudospin 
ideas following Ginocchio |TB|. In selecting Sp(3), we have been guided by experience with 



the noncompact Sp(3,R) algebra |^, the Interacting Boson Model and the Fermion 



Dynamical Symmetry Model [|T^]. In fact, our model contracts to the Interacting Boson 
Model as K oo. 

Results for the excitation energies of the lowest J = 0, 2, 4, 6, and 8 states are plotted 
in Fig. 1 for = 2k,, k = 5, 10, and 20, and coupling constants G = 0.05 and x = 0.005 
(dimensionless units) . The N/ k = 2 ratio corresponds to selecting the same one-third- filled 
shell nucleus in each case. The remarkable feature of the results is the sharpness of the 
transition from the superconducting to the rotational phase. For all k values shown, the 
system is essentially in a superconducting phase for a < 0.5 and a rotational phase for 
a > 0.6. The figure also shows the expectation values (S'+S'_)q, and {C2)a as functions of a 
for each J, where C2 is the SU(3) Casimir operator. The results exhibit the extraordinary 
purity of the wave functions on either side of the phase transition. For large k values, 

6 



{S+S^)a remains close to its a = value until a approaches a critical value oq. For a > ao, 
{S+S^)a falls rapidly to a small J-independent value. Similarly, {C2)a holds closely to its 
a = 1 value until a approaches ao from above after which it rapidly falls to a set of equally 
spaced values for each J. 

The purity of the wave functions is also indicated by the closeness of the results for the 
excitation energies to what one would get in first order perturbation theory starting from 
the appropriate limit. Thus, if one determines energy levels to first order in a, one obtains 

Ej{a) = Ej{0) 

+a[G{S+S.)o - 4x(C2)o + 3xJ{J + 1)] • (14) 

The linear expressions for the excitation energies, shown in the figure, indicate that the 
exactly computed results only diverge substantially from the first order expressions when a 
is relatively close to ao- On the other hand, if one determines energy levels to first order in 
(1 — a), one obtains 

Ej{a) = Ej{l) 

-(1 - a)[G'(5+5_)i - 4x(C2)i - SxAJ + 1)] . (15) 

The way in which the vibrational excitation energies drop precipitously just prior to 
the phase transition are reminiscent of RPA (random phase approximation) results for 



quadrupole vibrational excitation energies. As shown by Thouless |]T9[, within the framework 
of the RPA, the collapse of a vibrational energy to zero indicates an instability of the model 
ground state against vibrational fluctuations and the expectation of a phase transition to 
a deformed state. The present model, having a fermion foundation is clearly amenable to 
analyses along quasiparticle RPA lines. 

It is noteworthy that, for large k and = 2k, there is essentially no region of a for which 
the model is both rotational and superconducting. It appears that once the system starts to 
exhibit rotational states, the superconducting pair correlations are quickly suppressed. One 



recalls the similarity, noted by Mottelson and Valatin ||20|, between the Coriolis force acting 



on a particle in a rotating frame of reference and the interaction of a charged particle with 
a magnetic field. One knows that superconductivity tends to be excluded by a magnetic 
field; this is the Meissner effect. Mottelson and Valatin suggested that a similar collapse of 
superconducting pair correlations should be expected in rotational nuclei. However, their 
expectation was that the collapse would occur above a certain critical angular momentum. 
Thus, we would not have been surprised to obtain coexisting rotations and superconductivity 
at low angular momenta. Indeed, one of our motivations for considering the present model 
was to study superconducting flows in rotational nuclei. Thus, it is of considerable interest 
to understand what is special about our model that precludes this occurrence. [Subsequent 
calculations reveal that the strong suppression of superconductivity occurs in our model 
only when the number of nucleons N is set equal to 2k, as in the calculations reported here. 
We find that, for N > 2k., the pair correlations cause a mixing of SU(3) irreps in a highly 
coherent way such that the rotational character of the states is preserved. Results of these 
calculations will be reported in the more detailed analysisof the model to follow.] 

Outside of a narrow transition region, both phases of the system show a remarkable 
resilience to being perturbed. This we believe is due, in part, to the fact that the two phases 
are associated with incompatible dynamical symmetries. Put another way, the eigenstates 
of one phase are spread extremely thinly over the eigenstates of the other. However, consid- 
eration of other models indicates that this alone is not sufficient for a sharp phase transition 
to occur. For example, to see a sharp phase transition, it appears to be important to have 
many particles but only two-body interactions. 
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FIG. 1. The upper figures show the excitation energies for the Hamiltonian of eqn. |T]for different 
values of k and = 2k. The lower figures show the expectations {S+S-) and (C2) as functions of 
a for the corresponding k values. 
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